Lattice gas model in random medium and open boundaries: 
hydrodynamic and relELxation to the steady state. * 

Mustapha Mourragui ^ and Enza Orlandi ^ 

Abstract We consider a lattice gas interacting by tlie exclusion rule in the presence of a random field 
given by i.i.d. bounded random variables in a bounded domain in contact with particles reservoir at different 
densities. We show, in dimensions d> 3, that the rescaled empirical density field almost surely, with respect 
to the random field, converges to the unique weak solution of a non linear parabolic equation having the 
diffusion matrix determined by the statistical properties of the external random field and boundary conditions 
determined by the density of the reservoir. Further we show that the rescaled empirical density field, in the 
stationary regime, almost surely with respect to the random field, converges to the solution of the associated 
stationary transport equation. 

1 Introduction 

In the last years there has been several papers devoted in understanding macroscopic properties of non 
equilibrium systems. Typical examples are systems in contact with two thermostats at different temperature 
or with two reservoirs at different densities. A mathematical model of open systems is provided by stochastic 
models of interacting particles systems performing a local reversible dynamics (for example a reversible 
hopping dynamics) in a domain and some external mechanism of creation and annihilation of particles on 
the boundary of the domain, modeling the reservoirs, which makes the full process non reversible. The 
first question that one might ask for these systems is the derivation of the hydrodynamic behavior (law of 
large number) for the locally conserved field in the non stationary and stationary regime. There has been 
important classes of models, see for example [ELS1,2] , [DFIP], [KLO] in which it has been proved the 
law of large numbers for the empirical density in the stationary regime. Typical generic feature of these 
systems is that they exhibit long range correlation in their steady state. These long range correlations have 
been calculated from the microscopic dynamics only in very few cases; mainly in the case of the symmetric 
exclusion process, see [Spl], the asymmetric exclusion process, see [DEL], and in the weakly asymmetric 
exclusion process, see [DELO]. More recently breakthroughs were achieved analyzing the large deviations 
principle for the stationary measure. We refer to [BSGJL] for a review of works on the statistical mechanics 
of non equilibrium processes based on the analysis of large deviations properties of microscopic systems. 

In this paper we focus on the first step. We derive the macroscopic limit in the stationary and not 
stationary regime (hydrodynamic limit) for a particles system evolving according to local- conservative 
dynamics (Kawasaki) with hard core exclusion rule and with rates depending on a quenched random field 
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in a cylinder domain d > 3 in which the basis, denoted F, are kept at different densities. The restriction 
on the dimensions is only technical. We comment on this later. The rates are chosen so that the system 
satisfies a detailed balance condition with respect to a family of random Bernoulli measures (the random 
field Ising model at infinite temperature). To model the presence of the reservoirs, as in previous papers, we 
superimpose at the boundary, to the local-conservative dynamics, a jump dynamics (creation and destruction 
of particle). The rates of the birth and death process depend on the realizations of the random field and are 
chosen so that a random Bernoulli measure with a suitable choice of the chemical potential is reversible for 
it. This latter dynamic is of course not conservative and keeps the fixed value of the density on the boundary. 
There is a flow of density through the full system and the full dynamic is not reversible. The bulk dynamic 
models electron transport in doped crystals. In this case the exclusion rule is given by the Pauli principle 
and the presence of impurities in the crystals is the origin of the presence of quenched random field, see 
[KW]. The presence of the random field together with the exclusion rule makes the problem high not trivial. 
The transport properties of such systems in the case of periodic boundary condition on F has been studied 
by Faggionato and Martinelli, [FM]. They derived in d > 3, the hydrodynamic limit and gave a variational 
formula for the bulk diffusion, equivalent to the Green-Kubo formula. They proved that the bulk diffusion 
is a deterministic quantity depending on the statistical properties of the random field. Later, Quastel [Q] 
derived in all dimensions for the same model investigated by [FM] the hydrodynamic limit for the local 
empirical density. Applying the method proposed by Quastel, we could extend our results in all dimensions. 
Since our aim is to understand the role of the randomness in the non stationary and stationary state and 
not the role of dimensions in the bulk dynamics we state and prove our results in d > 3. Dynamical Large 
deviations for the same model and always with periodic boundary conditions have been derived in [MO] as 
special case of a more general system discussed there. The bulk dynamics is of the so-called nongradicnt 
type. Roughly speaking, the gradient condition says that the microscopic current is already the gradient 
of a function of the density field. Further it is not translation invariant, for a given disorder configuration. 
In order to prove the hydrodynamic behavior of the system, we follow the entropy method introduced by 
Guo, Papanicolaou and Varadhan [GPV]. It relies on an estimate of the entropy of the states of process with 
respect to a reference invariant state. By the general theory of Markov Processes the entropy of the state 
of a process with respect to an invariant state decreases in time. The main problem is that in the model 
considered the reference invariant state is not explicitly known. To overcome this difficulty we compute the 
entropy of the state of the process with respect to a product measure with slowly varying profile. Since this 
measure is not invariant, the entropy does not need to decrease and we need to estimate the rate at which it 
increases. This type of strategy has been used in previous papers dealing with the same type of problems, see 
[KLO] and [LMS] , which considered generalized exclusion process of non gradient type. The main difference 
with the previous mentioned papers is the presence of the randomness in the model considered here. This 
forces to take on the boundary a jump process depending on the external random field. Important step 
to derive the final results is then a convenient application of the ergodic theorem. Further we show that 
the empirical density field obeys a law of large numbers with respect to the stationary random measures 
(hydrostatic) . This is achieved proving that it is possible to derive the hydrodynamic for the evolution of the 
empirical measures starting from any initial particle configurations distributed according to the stationary 
measure, even though it is not possible to identify the profile. Then we exploit that the stationary solution of 



21/<iecem!)er/2008; 9:30 



2 



the parabohc nonhnear equation is unique and is a global attractor for the evolution. These two ingredients 
allow to conclude. Similar strategy for proving the hydrostatic is used in the paper in preparation by Farfan 
Vargas, Landim and Mourragui, [FLM]. 



2 The model and the main results 

2.1. The model 

Wc consider the d— dimensional lattice with sites x = (xi, . . . , Xd) and canonical basis £ = {ei, . . . , e^} 
and we assume in all the paper that d>3. We denote by A := [-1, 1] x r'^"\ where T''"^ is the {d - 1)- 
dimensional torus of diameter 1 and by T the boundary of A. 

Fix an integer AT > 1. Denote hy Aj^ = {-N,- ■ ■ , N} x T^"^ the cylinder in Z''^ of length 2A'' + 1 with 
basis the {d — l)-dimensional discrete torus and by Tn = {x € An \ xi = ±N} the boundary of Ajv- 
The elements of Ajv will be denoted by letters x,y,. . . and the elements of Ahy u,v, 

For a fixed A > 0, let S^) = \—A,A]^'^ be the set of disorder configurations on Z'^. On S^j we define 
a product, translation invariant probability measure IP. We denote by IE the expectation with respect to 
iP, and by a = {a{x), x G ^'^}, a(x) G [— A, a disorder configuration in S/j. A configuration a G T,o 
induces in a natural way a disorder configuration ajv on A^r, by identifying a cube centered at the origin 
of side 2N + 1 with An- By a slight abuse of notation whenever in the following we refer to a disorder 
configuration either on A^ or on we denote it by a. We denote by <Sjv = {0, 1}"^" and <S = {0, 1}^"^ 
the configuration spaces, both equipped with the product topology; elements of <Sjv or <S are denoted by jj, 
so that r]{x) = 1, resp 0, if the site x is occupied, resp empty, for the configuration rj. Given a G T,jj, we 
consider the random Hamiltonian : Sn IR, 

^"(^) = - E "(^)^(^)- (2-1) 

We denote by ^^'^ the grand canonical random Gibbs measure on Sn associated to the Hamiltonian (2.1) 
with chemical potential X G IR, i.e the random Bernoulli product measure 

^n\v)= n { eK.HAl + J - (2-2) 

When A = 0, we simply write /i^. We denote by and when A = 0, the measure (2.2) on the 

infinite product space S. Moreover, for a probability measure ji and a bounded function /, both defined 
on S or Sn, we denote by E^(/) the expectation of / with respect to /x. We need to introduce also the 
canonical measures Vp'^ , 

xeAiv 

for p e [0, •••)!]• It is well known [CM] that the canonical and the grand canonical measures are 
closely related if the chemical potential A is chosen canonical conjugate to the density p, in the sense that 



21 /december/ 2008; 9:30 



3 



the average density with respect to ^'^^ is equal to p. So as in [FM] one can define the random empirical 
chemical potential and the annealed chemical potential Xo{p)- To our aim it is enough to consider Ao(p). 
For p e [0, 1], the function Xo{p) is defined as the unique A so that 



IE 



= IE 



ga(0)+A 
1 + e«(0)+^ 



= P 



(2.3) 



We will consider as reference measure the random Bernoulli product measure ^'^(^ on Sn defined for positive 
profile p : A ^ (0, 1) by 



^p(-) M = 11 I g[a(cc)+Ao(p(x/JV))] ^ I \ ' 



(2.4) 



if p(-) = p is constant, we shall denote simply u'^^^ = v^'^ . We denote by ry^'^ the configuration obtained 
from rj by interchanging the values at x and y: 



' r]{x) if z = y 
r}{y) if z = X 
^ r}{z) otherwise, 



(2.5) 



and by if the configuration obtained from r] by fiipping the occupation number at site x: 
Further, for / : Sm IR, x,y £ Ajv, we denote 



r]{z) if z ^ X 

1 — ri{x) if z = X. 



(2.6) 



The disordered exclusion process on A^r with random reservoirs at its boundary Fjv is the Markov process 
on Sn whose generator Cn can be decomposed as 



where the generators £%, C]^ act on function / : Sn — * iR as 

{c-%f) {v) = Yl E ^^(^' ^ + ^) [(Vx,x+e/)(?7)] 

where e is a generic element of 5, the rate Cat is given by 

{X, y; V) ^C^x, y; n) = exp { - ^^^'^^"^^^^ } ; 



(2.7) 



(2.8) 



and 



in) = E C\x/N,^)[f{r^^) - /(,?)] . 
xeFjv 



(2.9) 



(2.10) 
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To define the rate C^{x/N, r/) we fix a function b{-) on F, representing tlie density of the reservoirs. We assume 
that b{-) is the restriction on F of a smooth function 7(-) defined on a neighborhood y of A, 7 : ^ — > (0, 1) 
and 7(w) = b{u) for u e F. The rate C** is chosen so that £^ is reversible with respect to u^^^ 

C {x/N,ri) = ?7(x)exp| ^ + (1 - 77(0;)) exp | . (2.11) 

The first term in (2.11) is the creation rate, the second one is the annihilation rate. Next we recall the 
relevant properties of C'N{x,y]'i]): 

a) detailed balance condition with respect to the m,easure (2.2), 

b) positivity and boundedness: there exists a > such that 

a-^ <CN{x,y;v)<(i, (2.12) 

c) translation covariant: 

C%{x,v\ri) =C^'°'{x- z,y- z;tzv) ^ t,,C%{x - z,y - z\ri) , (2.13) 

where for z in Z"^, denotes the space shift by z units on iS x S^i defined for all ry G iS, a G Sd and 
5 : 5 X ^ iR by 

{■TzV){x) =r]{x + z), {Tza){x) =a{x + z), (Tzg){rf,a) = g{Tzr],Tza) . (2.14) 

We omit to write in the notation the explicit dependence on the randomness a, unless there is an ambiguity. 
The process arising from the full generator (2.7) is then a superposition of a dynamics with a conservation 

law (the Kawasaki random dynamics) acting on the whole An and a birth and death process acting on F. 
Remark that if b{-) = 60 for some positive constant 60, then the generator Cjy, see (2.7), is sclf-adjoint in 
L^(i'^J^) and the measure I'Jj'j^ is the stationary measure for the full dynamics Cn- In the general case, 
when b{-) is not constant, since the Markov process on Sn with generator (2.7), is irreducible for all N >1, 
there exists always an unique invariant measure but in general cannot be written in an explicit form. 



2.2. The macroscopic equation 

The macroscopic evolution of the local particles density p is described by the quasi linear parabolic 
equation 

|p(0,-) = Po, (2-15) 

[p(t,-)|r = K) for i>0, 

where D{p) is the diffusion matrix given in (2.17), b{-) E C^{T) represents the interaction with the reservoirs 
appearing as boundary conditions to be imposed on the solution, sec its definition before (2.11), and po '■ 
A [0, 1] is the initial profile. The diffusion matrix is the one derived in [FM]. To define it, let * 

(S, = {g -.S X Ad ^ M; local and bounded } , (2.16) 

* A function g : S X Ad — > IR is local if the support of g, Ag, i.e. the smallest subset of such that g depends only 
on {{r]{x), a{x)) X e Ag}, is finite. The function g is bounded if SUp SUp„ \9{fl: <^)\ < ^O- 



21/december/200S; 9:30 



5 



and for g £ (K, Tg{r]) = X^^-g^d (txp) [r], a). The Tg{r]) is a formal expression, but the difference Vo,erg(?7) = 
rg(77°'^) — Tg{r]) for e e 5 is meaningful. For each p € (0,1), let D{p) = {Dij{p), 1 < i,j < d} be the 
symmetric matrix defined, for every a G M"^, by the variational formula 



1 

(a • D(p)a) = — ^ inf V IE 



,Ao(p) 



C°(0, er, 7?){aiVo,e,r/(0) + (Vo,eTg)(r?)}' 
where Ao(p) is defined in (2.3), x{p) is the static compressibility given by 



(2.17) 



X{p) = IE 



(2.16 



for a,b £ IR'^, (a • 6) is the scalar vector product of a and b and, recall, " " (•) is the expectation with 
respect to i_i°''^o(p) ^ ggg after (2.2), the random Bernoulli product measure on S with annealed chemical 
potential Ao(p). In Theorem 2.1 of [FM] it has been proved, for d > 3 and for p € (0, 1), the existence of 
the symmetric diffusion matrix defined in (2.17). Further it has been proved that the coefficients Di,j{-) are 
nonlinear continuous functions in the open interval (0, 1) and there exists a constant C > 1, depending on 
dimensions and bound on the random field, such that 



I 

C 



<D{p)<C^ PG(0,1) 



(2.19) 



where I is the dxd identity matrix. One expects the matrix !)(•) to be extended continuously to the closed 

interval [0, 1] and actually to be a smooth function of p, [KW]. Wc will assume all trough the paper that 
D{-) is well defined in [0, 1] and Lipschitz in the open interval. The diffusion matrix D{p) in a solid, in a 
regime of linear response, is linked to the mobility \ct{p)-, see [Sp], via the Einstein relation 



D{p) = ^<^{P)X{P)-'- 

The x{p) is a smooth function of p in [0, 1] and it can be easily proven from (2.18) that 



2P(i-p)<x(p)<p(i-p); 



-p{l- p)^<a{p)<CM^- P). 



(2.20) 



(2.21) 



where C is a constant that may change from one occurrence to the next. 



Weak solutions By weak solution of (2.15) we mean a function p(-, •) : [0, T] x A — > iR satisfying 
(IBl)peL2((o,r);ifi(A)) : 

j ds(^ j II Vp(s,m) f dw) < 00 ; (2.22) 

(IB2) For every function G{t,u) = Gt{u) in Cy{[0,T]x A ), where A=] - 1, lIxT''-^ and Cl'^{[0,T]x A ) 
is the space of functions from [0, T] x A to iR twice continuously differentiable in A with continuous time 

o 

derivative and having compact support in A we have 
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/ du{GT{u)p{T,u)-Go{u)p{0,u)}- [ ds ( du{dsGs){u)p{s,u) 
Ja Jo J a 

= - J ds^j duD{p{s,u))Vp{s,u) -VGsiu)^ ; 

(IB3) For any t e (0, T], Ti{p{t, •)) = b{-), a.e.. 
(IB4) p{0,u) = po{u). a.e. 

Notice that, since the original particle model cannot have more than one particle at a lattice site any 
solution p of (2.15) is bounded between and 1. The existence and uniqueness of the weak solution of (2.15) 
when (2.19) holds and D{-) is Lipschitz continuous for p G (0, 1), can be done using standard analysis tools. 
We refer to [LSU], chapter V or [DL]. Further, one immediately obtains by the characterization of H^^{A), 
see for example [E], page 283, that dtp G L'^{0,T;H~^{A)). Recall that H-^{A) is the dual of H^iA), i.e. 
the Banach space equipped with the norm 

= sup [{v, /) : ||/||ffi(A) < l} . (2.23) 

Stationary solution We denote by p the stationary solution of (2.15), i.e. a function from A — > [0, 1] so that 
p e H\A), for G G C2(A ) we have 

/ duD{p{u))Vp{u)-VG{u) =0, 

Ja (2.24) 
Tr(p(.)) = b{-), a.e. 



2.3. The main results 

For any T > 0, wc denote by {rit)te[Q.T] the Markov process on Sn with generator N'^Cn starting from 
r]Q = rj and by := its distribution when the initial configuration is rj. Wc remind that wc omit to write 
explicitly the dependence on a. The P^ is a probability measure on the path space D{[0,T],Sn), which we 
consider endowed with the Skorohod topology and the corresponding Borel a— algebra. Expectation with 
respect to is denoted by E^. If p^ is a probability measure on Sn we denote Phn{-) = J^^ P^(-)jU^((i77) 
and by E^n the expectation with respect to P^w. For t G [0,T],r] G Sn, let the empirical measure be 
defined by 

Trf (7?) = 7r^(dn;»7t) = ^ E %{x)S,/N{du) , (2.25) 

where (5„(-) is the Dirac measure on A concentrated on u. Since ri{x) G {0, 1}, relation (2.25) induces from 
P^A' a distribution Q^jv on the Skorohod space D{[0,T], A4i(A)), where A^i(A) is the set of positive Borel 
measures on A with total mass bounded by 1, endowed with the weak topology. Denote by A1J(A) the subset 
of A^i(A) of all absolutely continuous measures w.r.t. the Lebesgue measure with density bounded by 1: 

Ml{A) = {n € Mi{A) : n{du) = p{u)du and < p{u) < 1 a.e. }, 
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A4i{A) is a closed subset of A^i(A) endowed with the weak topology and D{[Q, T],A4i{A)) is a closed subset 
of Z)([0, T], A^i(A)) for the Skorohod topology. To state next theorem we need the following definition. 

Definition Given a Lebesgue absolutely continuous measure p{u)du € A4i{A), a sequence of probability 
measures {p^)n>o on Sn is said to correspond to the macroscopic profile p if, for any smooth function G 
and 5 > Q 



Theorem 2.1 Let d > 3 and assume that D[p) can be continuously extended to the closed interval [0, 1]. 
Let pb^ he a sequence of probability measures on Sm corresponding to the initial profile po. Then, LP a.s. the 
sequence of probability measures (Q^n)jv>o is tight and all its limit points Q* are concentrated on p{t,u)du, 
whose densities are weak solutions of the equation (2.15). Moreover if D{-) is Lipschitz continuous for 
p G (0,1), then (Q^n)jv>o converges weakly, as N '\ oo, to Q*. This limit point is concentrated on the 
unique weak solution of equation (2.15). 

Denote by v^'^ the unique invariant measure of the Markov process {rit)te[o,T] with generator N^jCn. We 
have the following: 

Theorem 2.2 Let d > 3, assume that D{p) can he continuously extended to the closed interval [0, 1] and 
Lipschitz continuous for p G (0, 1). For every continuous function G : A — > iR and every S > 0, 



with p{-) satisfying (2.24). 

3. Strategy of proof and basic estimates 

3.1. The steps to prove Theorem 2.1 

To prove the hydrodynamic behavior of the system we follow the entropy method introduced by [GPV] . 
As explained in Section 1, since the reference invariant state is not explicitly known, we compute the entropy 
of the state of the process with respect to a product measure with slowly varying profile 7(-). We prove 
in Lemma 3.8 that, provided 7(-) is smooth enough, G^ sufiices, and takes the prescribed value b{-) at the 
boundary, the rate to which the entropy increases is of the order of the volume, N"^, i.e the same order of 
the entropy and for finite time T this implies only a modification of the constant multiplying N'^. 

We divide the proof of the hydrodynamic behavior in three steps: tightness of the measures {Q^i^')n>i, 
energy estimates and identification of the support of Q* as weak solution of (2.15) with fixed boundary 
conditions. We then refer to [KL], Chapter IV, that presents arguments, by now standard, to deduce the 
hydrodynamic behavior of the empirical measures from the preceding results and the uniqueness of the 
weak solution of (2.15). We state without proving the first two steps, tightness of the measures and energy 
estimates. The proof of them can be easily derived from results already in the literature, which we refer to. 
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(2.26) 




(2.27) 



Proposition 3.1 (Tightness) For almost any disorder configuration a e So, the sequence ((^^^^n>i 
tight and all its limit points Q* are concentrated on absolutely continuous paths 7r(t, du) = p{t, u)du whose 
density p is positive and bounded above by 1 ; 

Q*|7r : 7r(t,du) =p(t,u)du| = 1 , Q*|7r : < p(t, u) < l| = 1 . (3.1) 



Tightness for non gradient systems in contact with reservoirs is proven in a way similar to the one for non 
gradient systems with periodic boundary conditions, see [KL], Chapter 7, Section 6. The main difference 
rehes on the fact that for systems in contact with reservoirs the invariant states are not product probabihty 
measures and some additional argument is required. This can be proven as in [LMS], Section 6. 

In the next step we prove that for almost any disorder configuration a € E73 , every limit point Q* of the 
sequence {Qfif^)N>i is concentrated on paths whose densities p satisfy (2.22). 

Proposition 3.2 For almost any disorder configuration a <E Sd, every limit points Q* of the sequence 

{Q^n)n>i is concentrated on the trajectories that satisfies (IBl). 

The proof can be done applying arguments as in Proposition A. 1.1. of [KLO]. However the latter proof 
requires an application of Feynman-Kac formula, for which we have to replace our dynamic (2.7) (cf. [FM]). 

We then show that IP— a.s. any limit point Q* is supported on densities p satisfying (2.15) in the weak 
sense. For (, e IN, x € Ajv, with —N + i<Xi<N — £ denote by r]^{x) the average density of 77 in a cube of 
width 2£ + 1 centered at x 

^ ^ v-\y-x\<e 

For a function G on A, e G £, d^G denotes the discrete (space) derivative in the direction e 

{d^G){xlN) = N[G{{x + e)/N)-G{xlN)] with x and x + eeA^r, (3.3) 
and to short notation we denote by G := d^^G for 1 < fc < d. 

Proposition 3.3 Assume that D{p) defined in (2.17) can be continuously extended in [0,1]. Then, for 
almost any disorder configuration a G Su, any function G in Cj['^([0,T]x A) and any 5 > Q, we have 

limsuplimsuplimsupPj^iv (jB^c^l > ^) = 0, (3.4) 

c— »0 a— >0 AT— »(x) 

where 

Blf = N-^ G{T,xlN)r,T{x)-N-^ ^ G{Q,x/N)r,o{x) - N-'' f dsG{s,xlN)r,s{x)ds 

xQKn xeAiv xeAjv ^ 

+ 



[ dsN'-'' Y i^^G){s,x/N){D,,^(r,^:^ix)) 



X {(2c)-^ [?7i"^l(a; + cNem) - V^^^^K^ ' cNem)]}} ■ 

(3.5) 
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The proof is given in Subsection 3.3. 

The last step states that IP— a.s., any Umit points Q* of the sequence ((5^iv)Ar>i is concentrated on the 
trajectories with fixed density at the boundary and equal to b{-): 

Proposition 3.4 IP— a.s., any limit point Q* of the sequence (Q|un)jv>i is concentrated on the trajectories 
that satisfy (IBS). 

The proof is given in Subsection 3.4. 
3.2. Basic estimates 

Lemma 3.5 (Ergodic lemma) Let V : J^d x A — > iR a bounded function, local with respect to the first 
variable and continuous with respect to the second variable, that is for any a G Sd the function u V{a, u) 
is continuous and there exists an integer i> 1 such that for all u & A the support of V{-, u) C {—i, ■ ■ ■ , i}'^- 
Then 

lim N-'^ V T^V{a,x/N)= [ IE\V{-,u)]du IP a.s.. 



(3.6) 



Proof. We decompose the left hand side of the limit (3.6) in two parts 

N-" r,V{a,x/N)=N-'' ^ (^T,V{a,x/N) - IE[V{-,x/N)]') 

+ N-'^ Y IE[V{-,x/N)] - I IE[V{-,u)\du . 

By the stationary of IP and the continuity oiu ^ IE\y{-, u)\ , the second term of the the right hand side of 
the last equality converges to as A'' — > oo. The first term converges to 0, from Chebychef inequality and 
the classical method of moments usually used in the proof of strong law of large numbers. □ 

We start recalling the definition of relative entropy, which is the main tool in the [GPV] approach. Let 
'^pI') product measure defined in (2.4) and /U a probability measure on Sn- Denote by H{n\u'^^^) the 

relative entropy of fi with respect to 

if(M|i^;(.^) = sup{y f{vMdv)-log j e^Wi.;('.^(d7?)} , 

where the supremum is carried over all bounded functions on Sn- Since t'pj-'.^ gives a positive probability to 
each configuration, n is absolutely continuous with respect to v'^'^^ and we have an explicit formula for the 
entropy: 

= f'og{-%,}d^. (3.7) 

Further, since there is at most one particle per site, there exists a constant C, that depends only on p(-), 
such that for all a 

H{n\u;i'^) < CN'^ (3.8) 
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for all probability measures /i on <Sjv (cf. comments following Remark V.5.6 in [KL]). 

It is well known that one of the main step in the derivation of hydrodynamic limit for the empirical 
density is a super exponential estimate which allows the replacement of local functions by functionals of 
the empirical density. One needs to estimate expression such as < Z, / in terms of Dirichlct form 

< —Cn\/ fir/), a/ f{ri) >^N, where Z is a local function and < •,• >^]v represents a scalar product with 
respect to some state . Since in the context of boundary driven process the invariant state is not explicitly 
known and we fix as reference measure some product measure v, see Lemma 3.6, there are no reasons for 

< —€,N\/f{ii), \/ f{rf) >v to be positive. Next lemma shows that this expression is almost positive. Let 
V%{-, v), V%{-, v) be functionals from h e L'^{v) to iR+: 



(3.9) 



xerjv 

Lemma 3.6 Let 7 : A ^ (0, 1) be a smooth function such that 7|p = For any a S S/j and a > 

there exists a positive constant Co = Co{A, ||V7||oo) so that for any f G L'^(^i'"l^) , 

f{v)C%f{r,)du^l'^{v) < -(1 - +C'o^'"'(«+ (3-10) 

/ f{v)C%f{v)du^l'^ iv) = • (3.11) 

JSn 



Proof. By (3.9) , 

/ f{v)^%f{v)du^l'; in) = -Kif, <(o ) 

+ E /Civ(x,x + e;7?)(V.,.+e/)(»?)/(»?"'"+^)^ii(a;,a; + e;7?)di.^('5^(7?), 

where 

i?i(x,x + e;r7) = (V.,.+er/(x)) (e^^-'^'^o^^^^/^))) - l) . 
By the elementary inequality 2uv < av? + a~^v'^ which holds for any a > 0, for any a;, a; + e e 

J Cn{x, x + e; v){V,,,+ef)f{v'''^+'')Ri{x, x + e, ri)du'^f^ (r?) 

-hj + + | / Cr,{x,x + e■,7l)f{7l^^-+'f{R^{x,x + e)fdu'^l'l{r^) . 

To conclude the proof it remains to use Taylor expansion and an integration by part in the second term of 
the right hand side of the last inequality. On the other hand, since 7]^ = 6(-) the measure v'^'^^ is reversible 
with respect to A simple computation shows that 

/ f{v)Cj,f{ri)dv;^^) iv) = -V%{f, ) . 
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□ 



Lemma 3.7 Let p,po ■ A ^ (0,1) be two smooth functions. There exists a positive constant C'q 
Cq{A, ||Vpo||oo5 ||Vp||oo) such that for any probability measure /j.^ on Sn and for any a G Hn, 



(3.12) 



Proof. Denote by f{'q) = and h{r]) = ^^{v)- Since f{r]) = h{ri) ^^"i.'v we obtain for e G f 



p(-) 



P0() 



p(-) 



and x,x + e G Ajv the following 

/ CN{x,x + e;ri) Vx,x+ey/f{v) dVpl^iv) 

J Sn 



I Cn{x, x + e;vi) [\/^(r?^'^+'=)ii2(a;, x + e;ri) + Vx,x+e\/^(??)] ^du^^^iv) 



<2 Cn{x,x + e;ri) Vx,x+e'^{n) 



Sn 



+ 2 / CN{x,x + e-,7^)h{ii-^-+%R2{x,x + e-^)Ydu;f^,^{r^), 

J Sn 



where 



R2{x,x + e;r,) = e^l,{{l/2)N-^d^[Xo{p{x/N)) - Mpo{xlN))]Wx,x+e'n{x)} - 1 • 
We conclude the proof using Taylor expansion and integration by parts. 



□ 



Denote by the semigroup associated to the generator N^Cn- Given a probability measures p on Sn 



denote by /i^(t) the state of the process at time t : ^i^{t) = ^^S^ . 

Recall that 7: A — > (0, 1) is a smooth profile equal to b at the boundary of A. Let be the density of 
IJ,^{t) with respect to ■ Let JO.*^^ be the adjoint of £jv in L'^{i^"l^). It is easy to check that 



dth^ 



(3.13) 



Notice that C* ^ is not a generator because v";^ is not an invariant mcasTirc for the Markov process with 



7(-) 



generator £jv. We denote by -ffjv(i) the entropy of iJ.^{t) with respect to ^'"(^! see (3.7), 



H^it) := H{„^{t)\u^l';). 



(3.14) 



Lemma 3.8 There exists positive constant C = C(||V7||oo) such that for any a > and for any a e T,d 

dtHn{t) < -2(l-a)A^2pO^(y^,^.";^)-27V2l?^(y^,^.";^) + -7V^ 



Proof. By (3.13) and the explicit formula for the entropy we have that 



dtHN {t) =N^ f H^Cn log (/if) du 

J Sn 



7(0 
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Using the basic inequality a( log b — log a) < — (^^/a — Vb) ^ + (6 — a) for positive a and b, we obtain 



dtH^it) < -2N^V%{Jh?,u'^J^)-2N^V],{Jh?,u'^J^) 



Since 7('u) = h{u) for ?x G F, v^'^ is reversible with respect to This implies that 



(3.15) 



•^ = 0. 



We shall now obtain a bound for C%h^dv'^'^ in terms of V^. Denote hy R : M ^ M the function 
defined by R{u) = e" — 1 — u . A standard computation shows that 



J Sn 

= ^'E E f CN{x,x + e-,r,)hi'{r,)R{N-'d^Xo{j{x/N))V,,,+,rj{x))d,.^l^^^ (3_^g) 
+ E {d^Ml{x/N))jw,,,+Mhnr})di^;'^){ri), 

where Wx,x+e{v) is the current over the bond {x,x + e) : 

Wj;,a:+e(r/) = Cn{x, x + e;r]) [r?(x) - r]{x + e)] . (3.17) 

We will often omit to write the dependence of Wx,x+e(j]) on ^ ^-iid i]. By Taylor expansion and the elementary 
inequality |-R(u)| < ^e'"', we obtain using the fact that 7 is smooth and is a probability density with 
respect to '^"(.^j that the first term of the right hand side of the (3.16) is bounded by C N'^ for some positive 
constant C. On the other hand integrating by part, applying the same computations as in Lemma 5.1 of 
[LMS], we obtain that there exists a constant Co = C(||V7||oo) so that for any a > 

jWx,x+ehfdv'^f^ < ^JcN{x,x + e;r])(Vx,.+e^^ d^;^^ + Co{a + N-'} 

for x,x + e G Ajv- □ 

For z e Ajv, M G IN denote by Am{z) the intersection of a cube centered at 2: e An of edge 2M+1 with 
Ajv, i.e 

AMiz):={z + AM}r\AN. (3.18) 

For probability measure on <Sjv, denote by the Dirichlet form corresponding to jumps in 

Am{z): 

T>%M^^'') = \ E JcN{x,x + e;rj)iVx,.+ef{v)fdu''{r,). (3.19) 
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Similarly, for 2: e Fjv define 'D\j ,y^) the Dirichlet form corresponding to creation and destruction of 
particles at sites in Fjv which are at distance less than M from z : 

T^mM^^"") = \ I C\x/N,rj){f{vn - f{il)fdu^{n) . (3.20) 

xerjvnAM(z) 

Fix any z G F^r denote by f^'^ the Radon-Nikodym derivative of /J.^ (t) with respect to v'^^^^j^y the random 
Bernoulli measure on Sm with constant parameter equal to Recall that we denoted by the Radon- 

Nikodym derivative of with respect to v'^'j^ and that h{j^) = "y{j^) for z G F. We have the following 

result. 

Lemma 3.9 Take M e IN , M < N. There exists a positive constant Cq = C(||V7||oo) depending only on 
7(-) such that for any z G Fjv 



The proof is similar to the proof of Lemma 3.7. 

3.3. Proof of Proposition 3.3 

We prove in this section Proposition 3.3. Let Q* be a limit point of the sequence {Qij,n)n>i and assume, 

without loss of generality, that IP— a.s., Q^n converges to Q*. Fix a function G in C^'^([0,T]x A). For 
a e flo consider the P^w martingales with respect to the natural filtration associated with {r}t)te[o,T]: 
Mf = Mf and AAf = A/f t G [0, T], defined by 

Mf = < Trf , Gt > - < TT^, Go > - /* ( < Trf , d,G, > +N^Cn < Trf , G« > ) rfs , 

Jo 



A/f = (Mf - ^* {iV2£^ ( < Trf , G, > - 2 < Trf , G« > N^Lj, < 7rf , G« >} . 



(3.21) 



A computation of the integral term of TV/' shows that the expectation of the quadratic variation of 
vanishes as iV t 0. Therefore, by Doob's inequality, for every 5 > 0, iP = 1, 



lim P^„ 



sup |Mf I >5 = . (3.22) 

0<f<T -I 



Thanks to (2.13) and since for any s G [0,T] the function Gg has compact support in A, a summation by 
parts permits to rewrite the integral term of as 

f < TT^,dsG, >ds + f [n^-^Y^ V (af G,)(x/iV)P^,,,+e.(%)|ds, (3.23) 
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where the current Wx,x+ek is defined in (3.17). To localize the dynamics define for any < r < 1 

-rN <xi< rN}, 



Ar = [-r, r] X r" \ ArN = {{xi, ■■■,Xd) e An 
TrJV = {x £ ArN '■ Xi = ±rN}. 

Set, for < a < c < 1, A; = 1, . . . , rf, 



(3.24) 



Vf''^'«(77,a) = NWo,e, + Yl (^'""^^0)) {(2c)"' [v^'^^'HcNem) - v^'^^'^-cNem)] } • (3.25) 

m— 1 

Next theorem is the main step in the proof of Proposition 3.3. 

Theorem 3.10 Assume that D{-) defined in (2.17) can he continuously extended in [0,1]. Then, IP = 1, 
for any G € Ci'2([0,T]x A), 



lim sup lim sup lim sup 

c->0 a^O JV->cxD 



N- 



J2 G,(x/7V)T,Vf'='^(r,„a)ds 







(3.26) 



for k = 1, . . . ,d. 

Proof. Lot < 9 < 1 such that for any t G [0, T], the support of the function Gt is a subset of A(i_29)- Fix 
a smooth fimction jg: A (0, 1) which coincides with b at the boundary of A and constant inside A(i_^). 
Denote by Z'^''"°'{G,r]) the quantity 

Since the entropy of fx^ with respect to i^^'f^ is bounded by C^IAjv] for some finite constant Cg, by the 
entropy inequality, the left hand side of (3.26) is bounded above by 



(3.27) 



for any positive B. Since e'^' < + e^^ and limsupiV^'^log{a7v + ^Af} < maxjlimsup A^^'^logaAr , 
limsup A^~''log6jv}, we may remove the absolute value in the second term of (3.27), provided our estimate 
remains in force if we replace G by —G. By the Feynman-Kac formula, 

where XN,c,a{Gs) is the largest eigenvalue of the TV^j/;^™ + BZf '^'"(G^, /?)} where r^"" := i(£jv + jO.*,^n) 
and -C*^ ^ is the adjoint of Cn in ^^(^7J(^))- -^y variational formula for the largest eigenvalue, for 
s e [0, T] , we have that 

-g^Ajv,c,a(Gs) =sup|y (G,,r?)/(7?)i/"^' (dT?) + — <CN\n,\l >ie(-) I • 
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In this formula the supremum is carried over all densities / with respect to ^'"j^!) and notice that we used 

< Cn\/J, -/J >7e(.)=< Cj^™'y/f, \ff >7e(.)- Since 7e(-) coincides with on F, C}^ is reversible with respect 
to 7e(-); so that < C^]^^T,Vy >je(-) is negative. Wc then apply simply (3.10) of Lemma 3.6 with a = 1 to 
estimate < Cj^^fJ ,\fj >7(,(.) by — (l/2)I?5v('\/7) ^"j^^)) + C'gN'^~^ for some constant C'g. In particular, to 
prove the theorem, we just need to show that 



limsuplimsuplimsup / dssup{ f ''='"(G,, 7?)/(r?)^"f Jdr;) - i^V^-dpO^^yj^ ^a,ArO ^ q 

for every B > and then let B '\ oo. Notice that for N large enough and a, c small enough, the function 
Z^''''°'{Ga,Ti) depends on the configuration rj only through the variables {?7(x), x e \i-e)N}- Since ^e{-) 
is constant, say equal to 70 in A(i_g), we may replace ^"'J^.y in the previous formula by f^^^ ■ The v^^^ is 
reversible for C% and therefore ©^(^ , ^'"o^) is the Dirichlet form associated to the generator C%. Since the 
Dirichlet form is convex, it remains to show that 

limsuplimsuplimsup dssup{ [ >^'"(G,, r,)/(r,)i."'|^.(dr?) - ^iV^-dpO^^^^ atO ^ q 

for every B >Q. This result has been proved in [FM], Theorem 3.2. □ 
Proof of Proposition 3.3: By (3.21), (3.23) and (3.25), applying Theorem 3.10 we obtain (3.4). □ 
3.4. Proof of Proposition 3.4 
For a > 0, u G A denote 

iJu) = Irr -, (u); (3.28) 

and for A c A define the sets A"^ as 

A+ = {{ui,...,Ud) e A : ui>0}, A- ={{uu...,Ud)GA : ui<0}. (3.29) 

Wc define similarly A+ and A]^ when An C An- Let G(-,-) G G^'^ ([0,T] x A), G D([0, T], 7Wi(A)) and 



1^ Cllil^ 

for < a < c < 1, define the following functional 

T 

is I 

A(l-c) 



-P^c(m(t)) = / d.s dulGsiu) {2c) ^ (fis * i-a)iu + cei) - (iis * i'a.){u - cei)\\ 

Jo Ja(1-c) ^ L /„ „„N 

^ ^ (3-30) 

+ j ds j dude-,Ga{u)[lJLa-k ba){u) - j ds| j b{u)ui{u)G a{u)d&^ , 

where Gs{u) = G{s, u), n=(ni, . . . , n^) is the outward unit normal vector to the boundary surface F and dS 
is the surface element of F. The proof of Proposition 3.4 follows from the next lemma. 

Lemma 3.11 For G(-, •) G G^'^ ([0,T] x A), IP a.s. we have 



lim sup lim sup lim sup 

c— >0 a— >0 JV— >oo 



<c(m^(-,-)) 



0. 
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Proof. To short notation, denote /«(«) := {Hs * La){u)- Taylor expanding we have that 
/ du\Gs{u) {2cy^ fs{u + cei)- fs{u-cei) \ 

JA(l-c) L J 

= 77-/ Ga{u- cei)fa{u)du- ^ [ Gs{u + cei)fa{u)du 

- / deiGs{u)fs{u)du + c / R{G,c,s,u)fs{u)du. 

JA(l-2c) Ja{1-2c) 

where \R{G,c,s,u)\ < sup^^A supg£[o^7'] \d^^Gs{-)\- Since < 1 uniformly in s and u 



(3.31) 



A(l-c) 



i?(G, c, s, u)fs{u)du 



<2sup sup \d^^Gs{u) 
ueAse[o,T] 



(3.32) 



and 



/ de^Gs{u)fs{u)du- / de-,Gs{u)fa{u)du 
Jk{1-2c) Ja 



< 2c sup sup |5eiGs(w)|. 
ueAse[o,T] 



Taking in account (3.31), (3.28) and (3.32) the lemma is then proven once we show that IP = 1 the following 
holds 



lim sup lim sup lim sup E^? 



a;e(A(i_„)jv\A(i_o_2c)N)* 



-]vi:rEH^)G4)}|]=o. 



(3.33) 



where for < £ < 1, A^at and (A^iv)"'" are defined in (3.24) and below (3.29). By adding and subtracting 
the same quantity in the expectation of (3.33), it is easy to see that the limit (3.33) follows once the next 
two lemmas are proven. □ 

Lemma 3.12 For G{-, •) e C^'^ ([0,T] x A), IP a.s. we have 



lim lim sup lim sup lim sup E w 

£^00 c^O a->0 N^oo 



Gs{x/N)rjt^{x) 



a;e(A(i_„),^-\A(i-a-2c)jv)=' 



-]V^ E Gs{x/N)vi{x)}\]=0. 



(3.34) 



Lemma 3.13 For G{-, •) e C^'^ ([0,T] x A), IP a.s. we have 



lim lim sup E A 

e^oo Ar->oo 



(l-T^)JV 



^ E &(^/A^)a(a=/A^)}|]=0. 



(3.35) 
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Proof of Lemma 3.12. The summation in (3.34) contains two similar terms. We consider the one 
corresponding to the summation of the right hand side of A;v (i-e. the one with signe +). By Taylor 
expansion applied to the function G, the expectation in the statement of the lemma is bounded above by 



E„ 



N{l-a) 



"'O Serl-i xi=iV(l-a-2c) + l 



where for xi G [—-/V, N], x — (x2, • • • , Xd) G ^ the vector (xi, x) stands for the element (a;i, X2-i ■ ■ ■ , Xd) € 
Ajv. We denoted by R{N, a, c, G) a quantity so that for G e C^'^ ([0, T] x A), 



limsuplimsuplimsup |-R(A'', a, c, G)| = 0. 

c— >0 a—yO JV— >oo 



(3.36) 



The next step consists in replacing the density average over a small macroscopic box of length aN by a large 
microscopic box. More precisely, for N large enough the expectation of the last quantity is bounded above 
by 



C||G||oo sup E^. 

2e<\y\<2Nc 



ds 



1 



T E k{{N-e,x)+y)-7ii'\N-e,x)] + R{N,a,c,e) , (3.37) 



where for all i, R{N, a, c, t) satisfy (3.36) and G is a positive constant. Observe that the first term of the 
previous formula is not depending on a but only on c, N and 

In view of the estimate (3.12) and Lemma 3.8 on the Dirichlet form and the entropy, by the usual 
two blocks estimate, the first term of (3.37) converges to an A'' | oo, c J, and £ t oo. That concludes the 
proof of Lemma 3.12. □ 

Proof of Lemma 3.13. The summation in (3.35) contains two similar terms, we consider the one 
corresponding to the summation of the right hand side of A;v It is easy to see that the expectation in (3.35) 
is bounded above by 



IIGIIc 



(3.38) 



For any fixed positive integer (. denote by = {(0, x) : x G T^j^'^, \x\ < £} = ({0} x T^"^) n A^(0), for 
notation see (3.18). For u gF, denote 

^lo(/>^) = ^ E/ C'oiu,x,r,){f{vn-m)'duir]), 



where 



Cl{u,x,ri) =77(a;)exp| - 



xer^ 



a{x) + \o{b{u)) 



} + (l-??(a;))exp| 



a{x)+Xo{b{u))- 



(3.39) 



The difference with the rate in (2.11) is that here u is fixed. Let ^^^) the product measure, see (2.4), 
where p{jj) = b{u) for Va; e Ajv and the restriction of to {0, Ij^^C^). Let / : <Sjv — * IR, denote by 
the conditional expectation of / with respect to the cr-algebra generated by {r]{z) : z G A^(0)} : 



,(z)=«.), .eAm}fi^)d,,^(u)iv) for all ^ G {0, 1}^^(°). 



21/<iecem!)er/2008; 9:30 



18 



Note that r/^(0) — b{u) depends only on coordinates on the box A^(0), then by Fubini's Theorem, 

ds\vi{y-ier)-b{y/N)\] =T j |r?^(0) - 6(y/iV)| (r_(,_to)/|''^)'Wc;<(i)(»?) (3-40) 

where /^'^ = y /^^ fs'^ds and for all < s < T, /J'^ is the density of with respect to the product mea- 
sure ^'^('^) with constant profile The density (T_(y_^ei)/T stands for the conditional expectation 
of T_(j,_£ei)/T' with respect to the cr-algebra generated by {r]{z) : z G A^(0)}. 

Remark that, since the Dirichlet form is convex and since the conditional expectation is an average, 

i\b,jT ( fi Iy,N\t \ ^ i\b,jT ( I jy^N a,N \ 

^e,0 [y {.T-(y-eeOfT ) > ^b(v/N) ) ^ ^e,0 [y T-{v-ee^)fT > ^b{v/N) ) 

- ^e,y-eei [y It > ^b{v/N) ) 

J 

Applying Lemma 3.9 we obtain from (3.41) 



(3.41) 



IT J '''b{y/N)J 



< 



< 



(3.42) 



for some constant Ct that depends on T. By the same argument we obtain the bound on the Dirichlet form 

d9 



eft' 



(3.43) 



For fixed and large enough, there exists a constant Ct, such that for all positive integer fc > 1, applying 
(3.42) and (3.43), we can bound by above the expectation (3.38) as following 



This last expression is bounded above by 



TllGIUATi-'^ J2 sup { / \rj\0) - b{y/N)\firj)di^^^^/^){rj) - k 2^?,o(v7, </;v)) 



(3.44) 
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where, for u e F, 



Further, since the function 



A^ = {f: />o, y/(Oc;<)(o = i}- 

is continuous on F, from Lemma 3.5, for all positive integers £ and k, the limit when A'' | oo of the expression 
(3.44) is equal to 



r||G|| 



dulE 



sup 



Since / r]l{0) — h(u) fdv'^^^^irO < Cb for some positive constant Cb that depends on ||6||oo, the integral 
on F in the last expression is bounded by 



[duIE\ sup { / r?'(0)-6(«)/(^X('iW} 



where for a positive constant C, kc '^^ following set of densities, 

C 



We first consider the limit when fc t oo and use the usual technics in the replacement lemma. Since for any 
£ > 1, any constant C > and any u gT the sets A^f. ^ are compacts for the weak topology, for all £ > 1 



lim sup sup 



v'{0)-b{u) mdu^(^){v) 



sup 



where 



Ale = {/ e A-, , ( v7, <i) = , V% (/f , ^«;^^) = o} . 



By dominated convergence theorem, it is then enough to show that. 



lim sup IE 



sup 



ri\Q)-b{u) fiv)du^(^^{v) 



. 



Now, it is easy to sec that, due to the presence of the jumps of particles in the Diriehlet form I?^q and the 
presence of the creation and destruction of particles in D^^ the set Af;^ = {!}. Thus, to conclude the proof 
of the lemma, it remains to apply the usual law of large numbers. □ 

Proof of Proposition 3.4 Let Q* be a limit point of the sequence (Q^iv)Ar>i and let {Q^Nk)k>i be a sub- 
sequence converging to Q* . By Lemma 3.2 Q* is concentrated on the trajectories that are in L^([0, T]; iJ^(A)). 
ForO<c< land for e Li([0, T], A1?(A)), such that /Lt(i, dw) = p{t,u)du with p{- , ■) e L^{[0,T]; H^{A)), 
denote by F^{ij,) the functional 

-^c^(m(->-)) = / du\Gs{u) {2c)~^ p{s,u + cei) - p{s,u- cei)\\ 

J ds j dude^Gs{u)p{s,u) — J ds^ J 6(u)ni(u)Gs(u)ds| . 



+ 
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Prom Lemma 3.11 and the continuity of the function fj, — > F^^{fi), we have 



lim sup 



sup J 

c^O 



Ff(/x) =0. (3.45) 



On the other hand, an integration by parts and Taylor expansion up to the second order of the function 
Gs{-) permit to rewrite as 

Pc{^^{■r))= I -t- I Gs{u)p{s,u)duds- [ ^ [ Gs{u)p{s,u)duds 

Jo J(A\A(i_2c))+ -^0 J(A\A(i_2^))- 



6(u)ni(w)Gs(w)dS + R{c) 



where i?(c) = R{G, c) is a function vanishing as c | 0. Further one has, see Theorem 5.3.2. of [EG], that 



lim I p. ^ A I / P{s,y)dy = Tr{p{s,u)) a.e uer,Vs, 

r^O \B{u,r) nA\ JB{u,r)nA 



(3.46) 



and then by dominated convergence theorem 

i-T 



lim (/.(.,•)) = ^ dsj^ {Tr{p{s,u)) - K«))ni(w)G«(«)dS 
This together with (3.45) implies 



T 

ds 

JT 



(TT{p{s,u))-b{u)^ni{u)Gs{u)dS 



which concludes the proof. □ 

4 Proof of Theorem 2.2 

The main problem in proving Theorem 2.2 is that we cannot associate to the stationary measure Vg"'^ a 
macroscopic profile according to definition (2.26). If this would be the case the result would be a corollary 
of Theorem 2.1. Denote by := Q^Cn the probability measure on the Skorohod space D(\0,T], Ai^ 

induced by the Markov process (tt^) = {TTNivt))^ when the initial measure is Vg"''^ ■ Denote by At C 
D(^[0,T],M) the class of profiles p{-,-) that satisfies conditions (IBl), (IB2) and (IBS). The first step to 
show Theorem 2.2 consists in proving that all limit points of the sequence ( Q^) are concentrated on At' 

Proposition 4.1 The sequence of probability measures ( Q^) is weakly relatively compact and all its 
converging subsequences converge to the some limit Q* that is concentrated on the absolutely continuous 
measures 'K{t,du) = p{t,u)du whose density p satisfying (IBl), (IB2) and (IB3). 

The proof of Proposition 4.1 follows the same steps needed to show Theorem 2.1. We just have to show 
the analogous of Lemmas 3.6, 3.7, 3.8 and 3.9 when the measure p,^ in the statements of these lemmas is 
replaced by Vs"'^ ■ The only lemma to be slightly modified is Lemma 3.8, see Lemma 4.2 given next. Recall 
that 7: A ^ (0, 1) is a smooth profile equal to b at the boundary of A. Let be the density of Ug"'^ with 
respect to the measure ■ 
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Lemma 4.2 There exists positive constant C = C(||V7||cxd) depending only on 7(-) such that for any a > 



C 



7(-) 



Proof. By the stationary of Vs°''^ , 

dtHr,{t) = f h^Cr, log = . 

Recalling that the generator has two pieces and applying the basic inequality a (log 6 — log a) < —[^/a — 
Vb)^ + [b — a) for positive a and b, we obtain 

0= / h^c^iog{h^)du;i!; < -2N^V%{Vh^,u^l!;) -2N^V>},{Vh^,u^l';) 



rO uN J a,N , ,t 
Lj^h +N 



2 ' c%h^d.;i''^ 



'Sm JSn 

We then apply the same computation as in the proof of Lemma 3.8, ((3.15) and (3.16)). 
Proof of Theorem 2.2 



□ 



Let Q* be a limit point of ( Q^) and ( Q^*") be a sub-sequence converging to Q*. Let p be the 
stationary solution of (2.15), see (2.24). We have by Proposition 4.1 the following: 

lim qN^{a) := lim Q^" { (tt^ ,G) - {p{u)du,G) ] 

fc— >oo fe— ►oo L J 

= Q* { I (p(T, ■),G)- {p{u)du, G) I t{AT}{p) } 

< IIGIloo q*{||p(t,.)-p(.)||,iMt}(p)}}. 

Denote by p°(-,-) (resp. p^{-,-)) the clement of At with initial condition /9"(0, •) = (resp. p^{0^-) = 1). 
From Lemma 5.8, each profile p(-,-) G At is such that for all t > Q, a|m G A : < p^{t,u) < p{t,u) < 

^{t,u) < l| = 1 and a|u e A : p^{t,u) < p{u) < p^(t,u)| = 1, where A is the Lebesgue measure on A. 



P 

Therefore 



hm q^.ia) < \\G\\^ \\p%T,-) - p\T,-)\\^ , IP = 1. 



To conclude the proof, it is enough to let T | oo and to apply Theorem 5.1. 



□ 
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5 Appendix 

In this section we show the global stability of the stationary solution of (2.15). 

Theorem 5.1 Global stability. Let D{-) be Lipschitz. Let p{t,po) be the solution of (2.15) with initial 
datum po, < po{u) < 1, u G A, and p the stationary solution of (2.15). We have 

lim / \p{t, u) - p{u)fdu = 

for all p> 1 . 



The proof of the theorem is based on an extensive use of monotone methods, see [S] . We were not able to find 
the precise reference, so we briefly sketch it for completeness. We need to introduce some extra notation. 
Let C^'^([0, T] X A) be the space of functions from [0, T]x A to IR twice continuously differentiable in A with 
continuous time derivative. Denote by 

g:= {Gg C^'^{[0,T] X A),G{t,u) = Gt{u) pointwise positive, G{t,u) = 0,Vw G r,Vt G [0,T]} . 

It is convenient to reformulate the notion of weak solution of (2.15) as following. A function •) : [0,T] x 
A — > [0,1] is a weak solution of the initial-boundary value problem (2.15) if p G L^(0, T; if^(A)) and for 
every G G Q 



/ du{GT{u)p{T,u) - Go{u)po{u)} - ds du {dsGs){u)p{s,u) 
J A Jo Ja 

= ds{ j^duAij{p{s,u))^Gs{u)- j^Aij{b{u))dnMs,u)ds} 



(5.1) 



where A,j{p) = D,j{p')dp' . A function p+(-,-) : [0,T] x A ^ M is a weak upper solution of the 
initial-boundary value problem (2.15) if /9+ G I/^(0,T; iJ^(A)) and for all G G 5 we have 

? / ^^{J^'^^ ^''^ ")) ~ ^''^ ^P'^^^' u))dn,G{s, u)dS^ 



- [ du{GT{u)p+{T,u) - Go{u)p+{u)} - I ds j du{dsGs){u)p+{s,u) <Q, 
J A Jo Ja 

Ti{p+{t,-))>b{-) on r 
, p+{0,u) > poiu) ugA 



(5.2) 



A weak lower solution p~(-, •) : [0, T] x A ^ i7? is deflned reversing the inequality in (5.2). 

By a solution of the stationary problem (2.15) we mean a function p G H^{A) so that for all G G C^(A), 

pointwise positive vanishing on F 



J2 [j^duA,,,{p{u)))^G{u) - ^ Ai^,{b{u))dn,G{u)ds] = 



(5.3) 
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As before we define upper and lower solutions of the stationary problem (5.3). A function p+ is an upper 
solution for the stationary problem (5.3) if p+ e H^{A) and for all G G C^(A), pointwise positive vanishing 
onr, 

V| / duAij{p+{u)))^G{u)- [ Aij{p^{u))dr,,G{u)ds}<0, 

Lo ^-^A ^^'^ ^ (5.4) 

Tr(p+) > 6 on T, 

A lower solution of the stationary problem (5.3) is defined reversing the inequality in (5.4). 

To apply the monotone method we first show the following comparison principle. 

Lemma 5.2 Let (resp. p^) he a lower solution (resp. upper solution) of (2.15), dtp^ G L^(0,T; H~^{A)), 
for i = 1,2. If there exists s > such that 

X{u G A : p^{s,u) < p'^{s,u)} = 1 , 

where A is the Lebesgue measure on A, then for allt> s 

\{ugA: p^{t,u) < p'^{t,u)) = 1. 

Proof Take s <t <T and 5 > 0. Denote by the function defined by 

Fs{a):= ^I{o<a<6} + (a-^/2)l{a>5}, a G IR. 

Let Ai := ^^(r) be the set 

As = \{t,u)G[Q,T]x A : Q < p^{t,u) - p^it^u) < 5^. 

By definition Tr(/9^ ^ P^) < a.e. and therefore Tr(i^^(/9^ ~ P^)) = 0- Since p^ ( p^ )is lower (upper) solution 
of (2.15), we have that 

dr-^ j Fs{p^{t,u)- p^{t,u)^ ^ j du Fs(^p\t,u) - p^{t,u)j - J du Fs (^p\s , u) ~ p^ (s , u)^ 

< -S-^ J'^dr J rfuV(p^-p2)-{£>(pi)Vp^-£>(p2)Vp2} 
= -S-^ f dT [ duVip" - p2) • D{p^)V{p^ - p^) 

J s J As 

- 5-^ I dT f duV{p^ - p") ■ {D{p^) - D{p')]Vp'' . 

Js J As 

(5.5) 

Since D{-) is strictly positive, see (2.19), the third line of (5.5) can be estimated by above 

-i^*dr duy{p'-p')-D{p^)V{p'-p')<-^ jj, du\W^P'-p')\\^- (5-6) 
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Further, by the Lipschitz property of D{-) we have on the set As, supi<j \Dij{p^) — Dij{p'^)\ < M\p^ — 
p^l < MS for some positive constant M. By Schartz inequahty, the last Une of (5.5) is bounded by 

6-^MA [ dr [ du\\V{p^ - p'')f + 5MA-' f dr [ du\\Vp''f (5.7) 

J s J As J s J As 

for every A > 0. By (5.5), (5.6), (5.7) and choosing A = M~^C~^ to cancel the term in (5.6) and the first 
term of (5.7) we have 

j duFs(^p\t,u)- p^{t,u)^ - j duFs(^p\s,u)- p^{s,ufj 
< SC-^M^ J dr J du\\VpY ■ 



Letting J, 0, we conclude the proof of the lemma because Fs{-) converges to the function F{a) = a]Io>o as 
510. □ 

By Lemma 5.2 we immediately obtain the following corollaries. 
Corollary 5.3 Let mo : A ^ [0, 1] be a measurable function. There is a unique weak solution p{t,mo) of 
the equation (2.15) with initial datum tuq. 

Corollary 5.4 Let mo be a lower stationary solution of (5.3). Let p{t,mo) be the solution of (5.1) with 
initial datum mo then p{t,u) > mo{u) a.e in {u,t). 

The proof is an immediate consequence of Lemma 5.2 with p^ := mo and p^ := p. When the initial datum 
of solution of (5.8) is an upper stationary solution we have: 

CoroUciry 5.5 Let mi be a upper stationary solution of (5.3). Let p{t,m\) be the solution of (5.1) with 
initial datum mi then p{t, u) < m\{u) for t G [0, T] and u G A. 

Next we show that when a lower (upper) stationary solution mo (mi) is taken as initial datum, the 
corresponding solution p{t,mo) {p{t,mi)) is monotone nondecreasing (nonincreasing) in time. 

Lemma 5.6 Under the assumptions of Corollary 5.4 p{t,mo) is a nondecreasing solution of (2.15) for all 
te%T]. 

Proof: Corollary 5.4 implies that p(s,mo) > mo for all s > 0, since mo lower solution. Let p{t; p{s,mo)) 
be the solution of (5.1) starting at time t = from p(s,mo). Then p{t; p{s,mo)) > p{t,mo) since the initial 
datum p{s, mo) > mo. But p{t; p{s, mo)) = p{t + s, mo) by uniqueness of weak solution then p{t + s, mo) > 
p{t,mo) > mo. □ 

Lemma 5.7 Under the assumptions of Corollary 5.5 p{t,m\) is a nonincreasing solution of (2.15) for 

ie [o,T]. 

The proof is similar to the one of Lemma 5.6. 

Lemma 5.8 Let mo be a lower solution and mi be an upper solution of (5.3), mo(-) < mi(-) a.e in A, we 
have 

mo < p{t;mo) < p{t;m\) < mi e (0,oo) 
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The proof is an immediate consequence of the previous results. □ 

Lemma 5.9 Under the assum,ption of Lemma 5.8 the solutions p{t: mo) and p{t; mi) exist for all t € [0, oo) 
and they converge in I'^(A) for p £ [l,oo) to limits and p*{-), both solutions of (5.3). Further 

P*{u) < p*{u) a.e. 



Proof: Since p{t;mo) is nondecreasing in t and p{t;mo) < rui for any t > 0, p{t;mo) converges almost 
everywhere in A as t — > oo and G L°°(A). By the monotone convergence theorem p{t;mo) — s- /0*(-) for 
p e [l,oo). Next we show that solves (5.3). Take as test function in (5.1) the following function 

0{t)F{u)] F{u) > 0; C> (3{t) > 5 > 0; P'it) > 0, (u, e A x iR+ 

/? e C^(ii+), F e C^(A) vanishing at the boundary. Then for all i > 0, see (5.1), we have 



(5.9) 



du{l3{f)F{u)p{t,u) - l3{Q)Fiu)po{u)] - i dsl3'{s) i duF{u)p{s,u) 

Jo J A 

= Yl fdsP{s){ J duAMs,u))^F{u)- J Aij{b{u))dn,F{u)ds]. 

Divide by t the left and right side of (5.9) and then let t — > oo. For the left side we have 

ijy" du{(3{t)F{u)p{t,u) - P{0)F{u)po{u)} - dsl3'{s) J duF{u)p{s,u)^ ^ 0. (5.10) 

By continuity of A{-) and since by assumption limg^oo P{s) = /3(oo) > 

lim / dsp{s){ f duAi,j{p{s,u))-^F{u)- ( A,,,{h{u))dn,F{u)ds] 

Jo '^Ja Jr ' 

= P{oo)J2{ J^duAij{p,{u))-^Fiu) - j^A,^,{h{u))dn,F{u)ds]. 



By (5.10) we then obtain 



(5.11) 



P{<x)Y,[ j^duAiMu))^F{u) - j^A,^j{h{u))dn,F{u)ds} = 0. 
Therefore is a solution of (5.3). The same can be argued for p* . □ 



The proof of Theorem 5.1 is a simple consequence of Lemma 5.9 and the unicity of the stationary solution 
p* = p^ of (2.15). 
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